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ABSTRACT 

The origin of the high inclination of Uranus' spin-axis (Uranus' obliquity) is one of the great unanswered 
questions about the Solar system. Giant planets are believed to form with nearly zero obliquity, and it has been 
shown that the present behaviour of Uranus' spin is essentially stable. Several attempts were made in order to 
solve this problem. Here we report numerical simulations showing that Uranus' axis can be tilted during the 
planetary migration, without the need of a giant impact, provided that the planet had an additional satellite 
and a temporary large inclination. This might have happened during the giant planet instability phase described 
in the Nice model. In our scenario, the satellite is ejected after the tilt by a close encounter at the end of the 
migration. This model can both explain Uranus' large obliquity and bring new constraints on the planet orbital 
evolution. 
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— k Introduction 
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^ | Today, Uranus' obliquity (97°) is essentially stable. This 



2. Numerical experiment 



^S'due to the regularity of its orbital motion and to the slow 
precession motion of Uranus' a xis compared to the sec ular 
!— -frequencies of the Solar system (jLaskar fc Robutelll993l ). It 
is 1 sometimes believed that a great collision with an Earth- 
rSiized protoplanet could be the reason of Uranus large obliq- 
uity. But this straightforward scenario hardly explains the 
O^presence o f regular satellites orbiti ng Uranus in its equato- 
rial plane (jKorvcanskv et aLlll990|) . However, the presence 
f^yf satellites around a planet can increase its precession rate 
pf by a large amou nt depending on their mass an d their 
^oijbital parameters (iTremainel 119911 : iGoldreichl 119661 : IWardl 
U1975I : iBoue fc Laskari 12*006^ For example, with a satellite of 
f-~mass m = O.OIMjj, where Mjj is the mass of Uranus, the 
• increase can reach a factor 1000 (Fig.[l|. The maximal effect 
. £^s' obtained for a satellite located at about 50 Uranian radii, 
k>vjiich is actually in the region wher e a satellite has been 
^V jbredic te d by some formation models (|Mosqueira fc Estrada! 
d£003allbl : [Estrada fc Mosaueiral[200rl . For comparison, the 
most distant regular satellite of Uranus is Oberon, its mass 
is 3.45 x 10~ 5 Mr/ , and its distance from Ur anus' barycenter is 
23 Uranian radii (jLaskar fc Jacobsonll987t) . The interactions 
between spin-axes and secular motions of the planets are also 
strenghtened when orbital inclinations are high. Such condi- 
tions could be met du ring the planetary m igration. Indeed, 
in the Nice scenario (jTsiganis et al.l [2005T ). Jupiter and/or 
Saturn should have undergone close encou nters with the ice 
giants to reach their present eccentricities (jMorbidelli et al.l 
2009). These close encounters can raise the inclinations. 
Moreover, the additional satellite can be ejected during one 
of these encounters. We therefore propose that Uranus had 
an additional satellite and that its spin-axis was tilted during 
the planetary migration. 



The construction of such a scenario for Uranus tilting can 
be described in four steps. 
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Fig. 1. — Uranus precession frequency in presence of 
a heavy satellite. Uranus effective precession constant as 
a function of the distance of an additional satellite of mass 
m = 0.01 Mu (a), m = 0.005 M v (b ), and m = 0.001 My 
(c), where Mjj is the mass of Uranus (|Boue fc Las kar 2006) . 
For this calculation, the semi-major axis and the eccentricity 
of Uranus are set to the current values, and the satellite is 
assumed to have a circular orbit. All inclinations as well as 
the obliquity are set to 0. 



First, we simulated the Nice model (jTsiganis et a l. 2005). 
We integrated 10 000 migrations of the giant planets over 10 
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Fig. 2. — Comparison between orbital evolution and 
obliquity increase. a,b, Example of orbital evolution of 
the giant planets during the planetary migration over 2 Myr. 
(a) semi-major axis, minimum, and maximum heliocentric 
distances, (b) Uranus inclination, c, Maximal tilt starting 
from zero obliquity for any orbital evolution with the same 
semi-major axis, eccentricity and inclination as Uranus in 
subfigures (a,b). In this calculation (see Fig. [5]), the effect of 
an additional satellite at 50 Uranian radii is implicitely taken 
into account in the precession constant. We considered three 
different masses for the satellite : 10~ 4 Mjj, 10 -3 Mjj, and 
0.01 Mjj, where Mjj is the mass of Uranus. For each mass, 
the satellite eccentricity is set to (lower boundary) and 0.5 
(upper boundary), d, Results of a numerical integration with 
a satellite of mass m = 0.01 Mjj (black curve). The satellite 
is ejected by a close encounter with Saturn at t = 0.38 Myr. 
Once the satellite is ejected, the obliquity remains constant 
(gray curve). In these plots, the obliquity is measured relative 
to the fixed plane orthogonal to the final total orbital angular 
momentum. 




Myr. For these simulations, 

8a . . . 
a = — cxp (— t/T) 

T 

de/dt = -e/(2r e ) 
di/dt = -i/{2n) 

with r = 2 Myr, and r e = n = t/10 as in dLee et al.ll2007t) . a, 
e and i are respectively the semi-major axis, the eccentricity 
and the inclination of the giant planets; t is the time and 5a is 
the difference between the initial semi-major axis of a planet 
and its current one. The initial semi-major axis of Jupiter 
was set to 5.45 AU. The initial semi-major axis of the other 
planets were obtained randomly with a uniform distribution. 
The initial semi-major axis of Saturn was varied in the range 
8.38-8.48 AU. The initial order of the ice giants was modified 
compare to the current Solar system : the initial semi-major 
axis of Neptune was varied in the range 9.9-12 AU and the 
initial semi-major axis of Uranus was varied in the range 13.4- 
17.1 AU. The terrestrial planets are not taken into account 
in this study. 

Then, out of the 5 142 simulations that survived without 
ejection or planet collision, we selected those where the planet 
final order is the same as in the Solar system. We obtained 
1 995 different integrations. As the tilt requires high inclina- 
tion, we kept only the simulations where Uranus' inclination 
increases beyond a given threshold. We set this threshold to 
17° which limits the number of simulations to 31. Among 
these simulations, we rejected those where the closest en- 
counter between Uranus and any other planet is closer than 
50 Uranian radii. With our criterion, we finally selected 17 
simulations. One of them is displayed in figure [3^, Ox Wc 
call it the reference simulation. 

In a third step, we computed the maximal effect of an ad- 
ditional satellite on Uranus' obliquity in any orbital evolution 
with the same semi-major axis, eccentricity and inclination 
as in the reference simulation, regardless the conjugated an- 
gles. For that, we used the expression of the effective precess- 
ion constant as a func tion of the satellite orbital parameters 
(|Boue fc L askar 2006). Then, we computed the maximal tilt 
given Uranus orbital evolution (jBoue et al.ll2009f ) . In the cal- 
culations, the satellite is at 50 Uranian radii in both circular 
and eccentric orbits. Figure [2£ shows the maximal obliquity 
that has been reached in these simulations. The evolution 
shows clearly that the tilt can only occur when the inclina- 
tion is high. In the present case, a satellite with m = 0.01A% 
is still necessary for the obliquity to reach 97°. 

Finally, we integrated the evolution of Uranus' spin-axis 
and the additional satellite in the 17 selected simulations. 
Calculations of the evolution of Uranus' spin-axis take into 
account the gravitational torques exerted by the Sun, by the 
additional satellite and by all the other giant planets. For 
each of the 17 planet migrations, we performed 100 integra- 
tions varying the initial semi-major axis of the satellite by a 
small amount (15 meters). The final obliquity distribution 
is given in figure |3l In 644 cases, the obliquity does not ex- 
ceed 10° because the satellite is ejected at the first encounter 
before the increase of the inclinations. But, if the satellite 
survive the first encounter, as in 62% of the cases, then the 
obliquity can reach large values. Among the integrations in 
which the satellite is ejected before the end of the migration, 
there is a final obliquity larger than 60° in 220 cases and an 
obliquity larger than 90° in 37 cases. 
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Fig. 3. — Distribution of Uranus final obliquity. This 
is the result of 1 700 integrations of Uranus spin-axis with an 
additional satellite : 100 per each of the 17 selected migra- 
tions. In black, the cases where the satellite is ejected; in 
gray, the cases where the satellite still orbits Uranus after 2 
Myr (the end of the integrations with a satellite) . The first 
bin has been troncated for a better visualization, its value is 
644 (black) + 8 (gray). Among the simulations with ejection 
of the satellite, there is a final obliquity larger than 60° (resp. 
90°) in 220 cases (resp. 37 cases). 



3. Dynamics of the tilt 

Here we explain the evolution of Uranus' spin-axis during 
the tilt presented in the figure [2j The smooth evolution of 
the obliquity during the tilt (Fig. 0b) suggests that it is due 
to a resonance. In the following, we show that the tilt ac- 
tually occurs during a 1:1 spin-orbit resonance between the 
precession of Uranus' axis and the regression of the node of 
its orbit. The obliquity e is measured relative to the invari- 
ant plane orthogonal to the total orbital angular momentum 
at the end of the simulation (10 Myr). Traditionally, it is 
defined relative to the orbital plane, but as the inclination 
rises to high values, it is preferable to use the invariant plane 
in order to avoid artificial evolution of the spin-axis. At the 
end of the reference simulation, Uranus' orbital inclination is 
very small (0.0024°) and the difference between the two defi- 
nitions is sufficiently small to be neglected. Let k) be a 
base frame such that the x-y plane coincides with this invari- 
ant plane. We note w, Uranus' spin-axis, and n, the normal 
of its orbit. The obliquity e is thus defined by cose = k ■ w. 
Let (j) a and (p v be the angles measured positively from the 
reference direction % to the projections of w and n into the 
x-y plane, respectively. The evolution of ip = <p a — (p y — ir 
is not steady but describes plateaus during the tilt (Fig. 
phases //, and IV +V). This confirms the 1:1 resonance be- 
tween the precession of Uranus' axis and the regression of the 
node of its orbit. 

In order to have a full understanding of the tilt, we now 
give the equations of motion that will allow us to describe the 
spatial evolution of the spin-axis displayed in the figure (Hi- 
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Fig. 4. — Details on Uranus tilting, a, Evolution of the 
resonant angle ip = 4> a — <j) v — 7r where 4> a and 4> u are angles 
measured positively from a reference direction to the projec- 
tion of Uranus' spin-axis w and Uranus' orbit pole n into 
the x-y plane, respectively, b, Evolution of Uranus' obliq- 
uity relative to the invariant plane, c, Evolution of Uranus' 
orbital inclination, d, Evolution of Uranus' spin-axis. The 
coordinates are x = sine cos ip, y = sine sin ip, and z = cose. 
The thin black circles in the x-y plane, that correspond to 
the thin black lines in the x-z plane and in the y-z plane, 
represent the locations where the obliquity is 45° (inner cir- 
cle), and 90° (outer circle). Uranus tilting is characterized 
by three resonant phases labelled //, and IV+V separated 
by non-resonant evolutions labelled /, III, and VI. 



Let a be Uranus' precession constant including the effect 
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of the satellites, and v the regression frequency of Uranus' or- 
bital node. In the frame rotating around k with the regres- 
sion frequency v, the Hamiltonian describing the spin-axis 
evolution reads 



H = ■ w) 2 - v(k ■ w). 



(1) 



In this frame k is constant, whereas n varies due to the evo- 
lution of the inclination. We have 




and w 



(2) 



with x = sine cos tp, y = sine sin tjj, and z = cos e. One can 
go from these variables to IWard fc Hamilto"nl (|2004j) ones by 
a rotation of angle i around the second axis. The equations 
of motion associated to the Hamiltonian (p} are given by 



which gives 



dw 
~dt 



or equivalently 



with 



— = V W H x w 

at 



-a(n ■ w)n x w — vk x w 



U X Z, 




(3) 



(4) 



(5) 



uj z = a cosi(z cos i — x sin i) + v, ^ 
io x = asini(z cosi — xsini). 

These equations are a combination of two rotations. The 
first one is a rotation around the third axis with the angular 
velocity uj z . The second is a rotation around the first axis 
with the angular velocity uj x . 

Now we describe the evolution of Uranus' spin-axis dis- 
played in the figured! 

During phase /, from t = to t — 289 kyr, the inclination 
i and the obliquity e are small, the three axes w, n, and k 
are almost aligned, there is no evolution. 

In phase II, the variable -0 remains close to — ir/2 (Fig. 
|3Ji), there is thus a 1 : 1 spin-orbit resonance. In that case, 
the angular velocity uj z in (J5J5]) is negligible and it remains 
only the rotation around the first axis 



(7) 



As the axis w describes an arc of a circle in the y-z plane 
(Fig. HJi), the third component z decreases, and thus the 
obliquity e increases (Fig. 

In phase III, the angle evolves (Fig. |4ji), hence the 
resonance is broken and the angular velocity lo z becomes im- 
portant. On the other side, the inclination decreases (Fig. 
3}:). Thus, the angular velocity lo x oc sini decreases too. 
Hence, the rotation in the x-y plane dominates (Fig. QJi) 
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Fig. 5. — Constraints on the migration and on the 
mass of the satellite, a, Lower boundary on the time 
required to tilt Uranus by 97° as a function of the inclina- 
tion o f its orbit for different precession constants (|Boue et al.l 
2009). b, Uranus' precession constant as a function of its 
semi-major axis for different ma sses of the satellite in cir- 
cular orbit ()Boue fc Laskar|[2006[ ). In case of eccentric orbit 
of the planet, the semi-major axis au should be replaced by 
a>u \/l — For example, if during the planetary migration, 
Uranus' inclination remains above 10° at least 1 Myr, then 
the precession constant should be larger or equal to 5"- yr _1 
to tilt Uranus (a) . Such a precession constant can be reached 
with a satellite of mass m = 0.001Af[/ if ajj^l — efj is less 
than 10 AU during the period of large inclination (b). 
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and the obliquity e is almost stationary (Fig. 0J>). 
During phase IV, the angle i/j is stable for a second time. 
The spin-axis is thus once again captured in resonance and 
the equations of motion are the equations (O . The axis thus 
describes an arc of a circle in the y-z plane (Fig. HJl). This 
time, it starts with y > and z < 1 (ip « n/2) and goes 
toward y = and z = 1. As z increases, the obliquity de- 
creases (Fig. 0}j). The orbital inclination i is similar in phase 
II, and Jy (Fig. 0J:), so is the angular velocity oj x . 

When the spin-axis crosses the x-z plane, the angle ip 
jumps from ir/2 to — n/2 and the system enters in phase 
V (Fig. 0k). The evolution is similar to the one of phase 
II. The spin-axis describe an arc of a circle in the y-z plane 
on the y < side (Fig. 0)1). The obliquity increases (Fig. 
0)d). However, as the inclination is higher than in phase II 
(18°, see Fig. 0J:), the angular velocity ui x is larger and the 
obliquity evolves faster (Fig. 0Jd). 

In phase VI, the resonance is once again broken (Fig. 0Ji), 
the angular velocity uj x becomes negligible with respect to 
u) z . The spin-axis describes an arc of a circle in the x-y plane 
(Fig. 0)1), and the obliquity remains constant (Fig. 0})). 

After the tilt, the satellite is ejected. The precession con- 
stant a decreases by a factor close to 1 000, which gives 
uj z ss v and oj x ss 0. The equations of motion §5§ thus show 
that the spin-axis precesses around the third axis at the an- 
gular velocity v. The obliquity remains constant as one can 
see in the figure [2j 
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4. Conclusion 



We have shown that the current obliquity of Uranus is 
compatible with planetary formation scenarios predicting 
small initial obliquities without the need of a large collision. 
Moreover, we confirm the necessity of the close encounters 
used in the Nice model to recover the present eccentricity of 
Jupiter and Saturn. Additionally, we solve the problem of 
the m issing satellite around Uranus (jMosaueira fc Estrada! 
2003a. b). Although stellite formation theories are not at the 
stage that they can constrain the final mass of the satellites, 
we acknoledge that the satellite we have introduced may be 
too massive. Nevertheless, a less massive satellite of only 
0. 001 Mu can still be sufficient to tilt Uranus if the plan- 
etary migration timescale is larger than the one used here 
(see figure [5]) . Several recent studies actually suggest such 
a longer migration timescale (iMurrav-Clav fc Chiang] 120051 : 
iBoue et "ail l2009t iLvkawka et al.l I2009T) . The parameters in- 
volved in the formation of the Solar system are still too nu- 
merous to be able to derive precise estimates of the probabil- 
ity of the present scenario to occur, but with this additional 
satellite, we are able to propose a scenario that fits with the 
present late migration scenario as given by the Nice model. 
Depending on the progress made in satellite formation theo- 
ries and possible future variations in the migration scenarios, 
the elements given here could provide strong additional con- 
straints on the migration timescales. 
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